ON THE VALIDITY OF CERTAIN 
CLASSICAL ASSUMPTIONS IN 
AIRCRAFT STABILITY 

C. R KELLY AND G R. RUPP 




I 1 ,.H brv lionl 

Muiterey, (Jaluomia 



ON THE VALIDITY OF CERTAIN 
CLASSICAL ASSUMPTIONS 
IN 

AIRCRAFT STABILITY 



t>y 



GEORGE R. KELLY 
Lieutenant 
U. S. Navy 



GEORGE R. RUPP 
Major 

U. S. Marine Corps 



UNIVERSITY OF MICHIGAN 
ANN ARBOR, MICHIGAN 



1 June 1955 



SUMMARY 



The purpose of this report was the investigation of the symmetric modes 
of motion of an airplane and of the concepts which have been used to obtain 
relatively simpler expressions for the period of these motions. The utility of 
the Reeves Electronic Analog Computer in this study was investigated. 

Specifically, the concepts concerning the symmetric motions were applied 
to the general body axes equations of motion in different ways . After the 
application of these simplifying assumptions, the modified equations were 
solved for the period of the motion concerned. Numerical data for a repre- 
sentative jet fighter type aircraft at 40,000 feet, Mach number .75* and at 
sea level, Mach number .25 was utilized in the solution by the different 
methods. These values for the period were compared to the true values as 
determined from the complete equations of motion. 

The possibility of increasing the accuracy of the period obtained in 
one of the methods was examined. 

A major objective of this report was the investigation of the utility 
of the analog computer in solving problems of this type . Mathematical treat- 
ment, while not difficult, is cumbersome and slow. The analog computer, on 
the other hand, after having been supplied with information, achieves a solu- 
tion almost instantaneously. It was desired to examine the effects of changed 
stability derivatives and coupling between the equations utilizing the computer. 
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The conclusions reached as a result of this investigation are: 

1. The concepts concerning the short period symmetric motion of an air- 
plane when applied to the general equations provide an acceptable 
approximation to the period of that motion. 

2. Concepts concerning the period of the phugoid mode, when applied, 
provide accuracies ranging from 30 percent low to 12 percent high, 
depending upon the method of application. It was concluded, there- 
fore, that use of the existing simplifications does not accurately 
describe the period of the true motion because essential couplings 
are neglected. 

3. The Reeves Electronic Analog Computer provided rapid accurate solu- 
tions to many sets of differential equations which had to be solved. 

The utility of the computer in investigations such as this in unques- 
tionable . Variation of stability derivatives or coupling between 
equations could be accomplished by changing one or several resistances 
and the new solution was immediately obtainable . 

This study was conducted at the University of Michigan, Ann Arbor, Michi- 
gan, during the period February to June 1955 under the supervision of Professors 
M. A. Brull and J. D. Schetzer of the Aeronautical Engineering Department. 



INTRODUCTION 



The advent of flexible high speed aircraft has somewhat complicated the 
problems which must be solved by the aerodynamic ist and the structures spec- 
ialist. With these complications has come an increasing tendency to write 
the equations of motion of the airplane in the body axes system of coordi- 
nates . 

The purpose of this report was the investigation of the symmetric modes 
of motion of a representative modern airplane in the body axes system and in 
particular, the investigation of the validity of the simplifying assumptions 
which have been used to obtain rapid approximate solutions for the period of 
these motions . 

The assumptions which have been made are based upon the observations of 
Lanchester who was the first to call the long period motion the phugoid. 
Lanchester stated that phugoid motion takes place at essentially constant 
angle of attack and further that the pitching motion is so slow that the 
aerodynamic moment may be considered zero, or nearly so. When these assump- 
tions are applied to the wind axes equations of motion and additional terms 
which do not materially affect the period are ignored, an expression is 
obtained for the period of the motion of .138 V. The same results have been 
obtained from considerations of kinetic and potential energy interchange? 

The current literature does not contain examples of the application of 
the phugoid simplifying assumptions in the body axes equations of motion. 
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The short period mode of motion is defined as a motion in pitch of very 
short duration which takes place as the angle of attack, which has been changed 
due to a disturbance, returns to the equilibrium value 0 The motion is so rapid 
that it is assumed that the forward speed does not change and consequently, 
the motion is independent of changes in perturbations in the forward speed. 
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SYMBOLS 



V P 

X 

z 

M 

u 

w 

q 



\ 

x u 



Xw 



Zw 

Z q 

Mu 

Mw 



aircraft velocity ft/sec 

aerodynamic force in X body axis direction lbs 

aerodynamic force in Z body axis direction lbs 

aerodynamic moment about the Y body axis ft lbs * 

perturbation velocity in X body axis direction ft/sec 

perturbation velocity in Z body axis direction ft/sec 

perturbation velocity in Y body axis rad/sec 

dot over a symbol - derivative of symbol with respect to time 

operator for derivative with respect to time, — 

dt 

dimensional derivative of the aerodynamic force in the X body 

Sx 

direction with respect to the perturbation velocity u: 

ou 

lb sec/ft 
ow 

^ lb sec/ft 

^ lb sec /ft 
dw 

ib sec/ft 

— ft lb sec /ft 

du 

ft lb sec/ft 

Sw 

M ft lb sec/ft 

^ ft lb sec/radian 
dq 



axis 

sec/ft 
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angle of attack of M.A.C. relative to wind - degrees 
angle of pitch of X axis measured from the horizon - degrees 
steady state lihear velocity along X axis - ft/sec 
steady state linear velocity along Z axis - ft/sec 
mass in slugs 

moment of inertia about Y body axis - slug ft 
Mach number 
Lift coefficient 



Sum of forces in X body axis direction 



SYMMETRIC EQUATIONS OF MOTION 
(BODY - AXES SYSTEM) 



The six equations of motion break down into two sets of three equations . 
One set describes asymmetric motion. The other set describes symmetric motion. 
The linearized equations for symmetric motion, neglecting Xq, in body axes 
coordinates are: 

F x = ° = [X u -mX]u + [X w ]w + [-mW ^ , c ° s 9 1 ] q 
F z = 0 = [Zju + [Z^mXjw + [Zq+mUx-M-Sin^jq 
My = 0 = [MuJu + [Mw-M^^w + [M q - Iyy X ]q 



The determinant of the coefficients is set equal to zero and the stabi- 
lity quartic is obtained 

AX 4 + BX 3 + CX 2 + DX + E = 0 



where 



A = 1 
B = - 



Xu M q . 

— “ — - - — — Za+riJi) 

m m I yy mI yy q 



+ — — g sin ©! 



% 



Mq W x 

= (Z ^ XU) + ^ + ' Iyy 



mI yy 

Zq+mUi i 

+ m 2 Iyy + m 2 



Mq s i n © t w 

D = (Zu X w-X u Zw) - ^ X {K^ 1 -K a Z v ) + 

* " yy yy 



m^Iyy 

(Zq+mUi) X w Mu g cos 0! „ „ s 

- ^ + (M ^«) 

E = S mIyy " 1 ^ Z u- M u Z w) + ^Sl^ -11 (Xw^’^Xu) 



The equations of motion are developed in Ref. 1. 
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When the analog computer is to be utilized, it is convenient to rewrite 



the equations of motion as follows: 



ml = X u u + X w w - m W x q-Mg cos 0 X 

m& = Zjju + Z^w + [Zq+mUxJq - mg sin © x 

# 

Iyyq = MyU + MtfW + M#w + Mqq 





The computer circuit for the solution of these equations is included as 



Fig. 1. 


The airplane 


numerical 


data supplied to the circuit is controlled 


by positioning the potentiometers . Each potentiometer 


controls a specific 


airplane parameter . 


These are: 








i. 


m 


6. 


*w 


11. 


Mu 


2. 


\ 


7 . 


mW 2 . 


12. 


Z 

u 


3. 


m 


8. 


mg sin 


13. 


Zq ■+ mU x 


k. 


Z w 


9. 




Ik. 


Mq 


5. 


I yy 


10. 




15- 


mg cos 0! 
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APPLICATION OF CLASSICAL ASSUMPTIONS 



FOR SIMPLIFICATION OF THE SOLUTIONS 
TO THE MODES OF MOTION 

The roots of the stability quartic yield the characteristics of the 
short period and phugoid motions. Various simplifying assumptions con- 
cerning these motions have been introduced in the past. Accuracy obtain- 
able after the assumptions have been applied varies . Normally, the short 
period motion is so heavily damped and so rapid that it is of little con- 
sequence . It does become important if the period is within the normal 
human response lag and not heavily damped, primarily because a pilot may 
inadvertently reinforce the motion when he attempts to correct it . The 
short period mode is important when platform stability is a consideration. 
The phugoid motion is a long period slow oscillation with weak damping and 
of such little consequence that it may be negatively damped and not render 
an airplane unacceptable. Because of the apparent unimportance of these 
modes in piloted subsonic aircraft, approximations to their period and 
damping are justified and desirable . 

However, the importance of these modes of motion in pilotless missiles 
at very high speeds has not been completely investigated. When a relatively 
long slender missile encounters severe turbulence or a sharp gust, the pos- 
sibilities seem to be present for considerable control difficulties. These 
motions could intensify the difficulties . 

It follows then that it is important to have a method for the rapid, 
accurate determination of the symmetric modes . 
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SHORT PERIOD METHOD 



The short period motion is obtained by assuming that the motion proceeds 
at constant forward velocity (u = 0). This follows if it is reasoned that by 
the time the forward velocity has changed, the motion has damped out and has 
ceased to exist. The assumption permits a reduction of the equations of 
motion to essentially a coupling between the perturbations w and q. The 
equations of motion become: 

F z = 0 = [Z^mXjw + [Zq+mUx- s * n .. 9 i jq 
My = 0 = [M w +M^\]w + [Mq - Iyy \]q 



Solution of these equations as written involves the solution of a cubic . 
The equation becomes a quadratic if the perturbation gravity force is assumed 
very small and ignored. The quadratic is: 

AX 2 + BK + C = 0 

where : 



A = 
B = 



1 

Tvv HI 



Mw 

mI yy 



( Zq+mU x ) 



C 



£y/Mg 

mI yy 



^w 

mI yy 



( Zq+mU i) 



The utility of the analog computer is demonstrated by the fact that the 
computer circuit for the short period method is achieved by uncoupling poten- 
tiometers 11 and 12 in Fig. 1. This uncouples and removes from the circuit 
all contribution to the short period of the u perturbation velocity and the 
solution is obtained using the reduced Z force and moment equations. 

Ref. 2 shows a simplified short period method in the wind axes system. 
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PHUGOID METHOD 



The long period phugoid motion is usually obtained by applying Lan- 
chester's original assumptions. Ref. 3 states that in a phugoid motion the 
angular pitching movement of the aircraft is at all times very slow and this 
implies that the aerodynamic pitching moment is nearly zero. Hence, the 
angle of attack is at all times very nearly that corresponding to zero 
pitching moment. This implies nearly constant angle of attack. Ref. 3 
further assumes that a and is strictly constant . 

This investigation revealed that a choice is available for the method 
of enforcing zero change in angle of attack in the body axes equations. 

3y definition: 

Wt 

a = arctan — L 
Ui 



if dW x : 
and dU x - = 
then 



da da 
dQ = d^+d^ 

perturbation velocity in Z direction = w 
perturbation velocity in X direction = u 



. da da 
Aa = w + bid u 

dWx dU x 



Aa = 



w cos a - u sin a 



V- 



P 



PHUGOID METHOD I 

The expression for the perturbation angle of attack is usually 
simplified by assuming the equilibrium angle of attack is small, and 
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because of this assumption, that 

By this method, the perturbation a is related directly to the pertur- 
bation velocity w. 

When Lanchester’s assumptions are applied to the equations of motion 
based on this near equality for A a, and Aa = w = 0, the stability quartic 
reduces to a quadratic which gives an approximate solution* The quadratic 
is obtained by ignoring the moment equation and deleting all terms in w in 
the force equations. Two equations remain in the variables u and q. These 
are: 

F x = 0 = [X^mXju - [mW ^ 11 ^ C ° S ®l ]q 

F z = 0 = [Zj u + [Zq+mUx-rcg. ®J-]q 



The determinant of the coefficients set equal to zero yields: 

AX 2 + BX + C = 0 



Where: 

A = [Zq+mU x ] 

r z n 

B = "t x u (-£ + UJ + mg sin © x + Z^J 
C = g [X u sin 0 X - Z^ cos 0 X ] 

When the analog computer is utilized, the two equations are rewritten 
as follows: 

mu = X u il - mW x q - mg cos 0 x q 
(Zq+mU^q = mg sin © x q - Z u u 

The computer circuit is shown in Fig. 2. Specific airplane parameters 
are controlled by the potentiometers as follows: 

1. Zq + mU x 2. mg sin © x 
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5. 


m 


6. 




7. 


*u 


8. 


mg cos ©x 


9. 


mWi 







Switch in Fig. 2 is open. 

This method of simplification produces a period which is greater than 
the true period. 



PHUGOID METHOD II 



A second method for the reduction for the complete equations to obtain 
the period of the phugoid mode was developed by the writers during this in- 
vestigation. This method involves the application of the identical simpli- 
fying assumptions of method I, but the complete equality 

_ w cos a - u sin a _ q 
V P 

was used to enforce no change in angle of attack. Then w = u tan a was 
substituted in the equations and the third equation was ignored. The three 
equations reduce to two equations which are: 



F = 0 = [X u -m\ + X w tan Ct]u - [m W x +- 5& ^ ® l]q 
F z = 0 = [Zu+Z^tana-mXtanCtJu + [Zq+mU 1 - rri ^' s ^^ l]q 
The determinant of the coefficients set equal to zero yields: 

AX 2 + B\ + C = 0 



where : 

A = m [Uj+V?! tan OC + ?_£L] 

x, m 

B = - [(Zq+mU 1 )(-jj-+ ^51 tan a) + mg (sin ©i - cos ©x tan a) 
+ w i (Zu +Z w tan a )J 
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C = g [sin ©i (x w tan a + X u ) - cos 0 X (Zy^-Zy tan a)] 
The equations are rewritten as follows when the analog computer is used: 
mi = (X u +X w tan oc) li - m Wj/j - mg cos ©i q 
(Zq+m Ui)q = mg sin ©j. q - (Zu+Z w tan cc) u + mu tan a 
The computer circuit is shown in Fig. 2. Specific airplane parameters 
are controlled by the potentiometers as follows: 



1. 


Zq + mUi 


6. 


Zy + Zw tan cc 


2. 


mg sin ©! 


7. 


X u + X w tan CC 


3- 


m 


8. 


mg cos ©! 


4. 


m tan OC 


9- 


mW x 


5* 


m 







Switch in Fig. 2 is closed. 

This method of simplification produces a period which is less than the 
true period. 

Solution by Method II provides a period which is exactly equal to .138 

v 

V^, arrived at in wind axes in Ref. 2, and is exactly equal to 4.44 de- 
y g 

rived in Ref. 3* Thus Method II is the true application of Lanchester's 
Basic Assumptions in body axes . These assumptions are based on the pure 
interchange of kinetic and potential energy . 

Ref. 4 states, "To secure a lift proportional to the square of the speed, 
the wings must exert a constant lift coefficient, so that their incidence must 
be constant . This implies that the aeroplane is statically stable as regards 
pitching, that it has a negligible moment of inertia about the pitching axis, 
and the length between wings and tail is small compared to the radii of cur- 
vature of the flight path . " 
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Method II does not substantiate this statement since substitution of 
w = u tan a, (for the condition Aa = 0) in the moment equation causes My 
and My to be grouped in the same coefficient and these terms exactly cancel 
each other because they become equal in magnitude and opposite in sign. 

( |Mu | = |My tan Ct | ) . This leaves the moment equation without a static sta- 
bility contribution and produces the same effect as zero static stability. 

In fact, solutions of pairs of the simplified equations with the moment 
equation included as one of the pair show the motion to be aperiodic . This 
is the motion to be expected from an airplane with zero static stability. 

The condition of zero static stability may be solved quickly on the analog 
computer by positioning the two potentiometers which control Mu and My on 
the total airplane circuit to zero. Computer solutions for this case show 
aperiodic motion. 

It is apparent that the Lanchester idealization enforces neglect of the 
moment equation because as far as the phugoid mode is concerned, the moment 
equation becomes inconsistent with the two force equations. This requires 
that the two force equations be used to obtain the solution for the phugoid 
mode and that an error in period must be accepted. 

MODIFIED PHUGOID METHOD II 

It should be noted that Aa == 2L = 0, if applied in the body axes 

V P 

equations, produces a period which is larger than the true period; whereas. 
Method II produces a period which is smaller than the true period. 

Comparison of the determinantal coefficients resulting from both Aa 
approximations was made and the terms in each of the coefficients which 
largely effected the period were examined. It was noted that the gravity 
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contribution in the coefficient multiplying X to the first power (Method II) 
was completely destroyed by tan a. 

[mg (sin ©i - cos ©x tan a)] a = 0 X 

The classical explanation of phugoid motion is based on the interchange 
of kinetic and potential energy. The complete removal of all effects of 
gravity, therefore, seems illogical. In short, if the factor tan OC is re- 
moved from the term mX tan a (second equation), the gravity effects are 
restored and the solution for the period of the phugoid is accurate to with- 
in one percent. The resultant period is nearly the actual value. 

This result is an experimental result arrived at on physical grounds 
alone and perhaps may not be applied with such accuracy in general because 
it cannot be mathematically justified. As shown above, the idealizations 
of Lanchester as generally applied are also not mathematically justified. 
(When Aa = 0, the airplane is not simultaneously statically stable in the 

equations of motion.) If the additional reasoning concerning the gravity 
term is accepted as described on physical grounds, the coefficients of the 
quadratic describing the phugoid mode are: 

A = m [Ui+Wi + 

B = - [(Zq+mU 1 )(^.+^l tan Oi) + mg (sin ©x - cos ©x) 

+ W x (Zu+Z^ tan a)] 

C = g [sin ©x (X u +X w tan a) - cos ©x (Zq+Zv/ tan a)] 

The computer solution for Modified Method II is obtained by ■chaining the 
setting of potentiometer 4 in Fig. 2. (switch -closed) to remove the factor tan a. 

Numerical examples for Method I, Method II, and Modified Method II are 
included in this report. Method I is 12 percent high. Method II is seven 
percent low and Modified Method II is within one percent, all referred to 
the actual period. 
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PHUGOID METHOD III 



One other method for an approximate solution for the phugoid mode has 
been used. This method makes use of the fact that phugoid motion is essen- 
tially a coupling between u andjq. The method of simplification used is to 
delete the w coupling terms in the force and moment equations thus uncoupling 
the Z force equation from the system. The X force equation and the moment 
equation thus reduced are solved simultaneously. The equations are: 



F x =' 0 = [X u -mX]u - [mW^ 11 ^ c ° s 0 l ]q 
My = 0 = [ Mu ]u + [ M q - Iyy X ]q 
The resulting cubic equation is: 

AX 3 + DX 2 + CX +D = 0 



where : 



A = 1 



B = 

C = 
D = 



m 



Mq + lyy ^u 



!_ m Xyy -> 

[X u Mq + m W x Mu] 

L ®iyy J 

[ Mu g cos Q x 1 

L lyy J 



This method produces an error in the period of 25 to JO percent . 

Solutions on the Reeves Electronic Analog Computer show this error . 

The computer solution is obtained by deleting the connections to potentio- 
meters 6, 9, and 10 in Fig. 1. This effectively uncouples the Z force equa- 
tion. The most powerful term in reducing the period 

was determined to be [Mw]w. This was shown by step by step replacement of 
the deleted terms. When [MwJw alone was reconnected, the period was restored 
to the correct value . 



15 



One other idealization was applied to the equations of motion in an 
attempt to find a simple accurate representation of the phugoid mode. Because 
the phugoid is described as a long slow oscillation about a datum line, it 
seemed reasonable to assume that motions of this type would be characterized 
by very small accelerations in the direction normal to the flight path. This 
may be expressed as: 

w - U x q = 0 

This equation was solved for q and substitution for q was made in the equations 
of motion. Pairs of the equation were solved with varying results. Solution 
of the two force equations produced aperiodic motion; solution of the Z force 
equation and the moment equation produced a short period mode frequency; and 
solution of the X force equation and the moment equation produced a cubic . 

This cubic was not solved but known accurate values of the phugoid roots were 
substituted to ascertain whether or not the correct phugoid mode was contained 
in the cubic . It was found that the roots did not satisfy the equation even 
remotely and the method was discarded. 

It seemed logical to assume further that if the normal accelerations are 
small or zero, that the condition of zero angle of attack change would occur 
simultaneously. When this additional condition was imposed, the equations of 
motion reduced to three equations in u but all the solutions indicated aperio- 
dic motion. 



NUMERICAL RESULTS AND DISCUSSION 

The following examples are selected to facilitate comparison of the 
results obtainable by the several methods discussed above . The airplane 
upon which this is based is a representative jet fighter type aircraft. 

Two case are presented: first, the case of Mach number 0.75 at 40,00 feet; 
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second, the case of Mach number 0.25 a b sea level. Normal mathematical 
solutions are presented and then are verified by Reeves Analog Computer 
solutions for the first case. The second case will be presented as analog 
computer solutions only. 



CASE I 

AIRPLANE DATA AND VELOCITIES 



S = 600 sq ft 
m = ll80 slugs 
Iyy = 80,000 slug ft 2 
h = 40,000 ft altitude 
a = 2.48 degrees = ©! 
V p = 728 ft/sec 
Ua. = 728 ft /sec 
W x = 31.5 ft/sec 
M = 0.75 



X u = -5.19 lb sec/ft 

Xy = 48.8 lb sec/ft 

Z u = -70.8 lb sec/ft 

Zy = -807 lb sec/ft 

My- = -637 lb sec 

M^j = 10.08 lb sec 2 

M u = 27.5 lb sec 

Zq =-4.94 x 10 3 lb sec 

Mq = -128.4 x 103 ib ft sec 



Using this data, the stability quartic coefficients are; 



A = 1 
B = 2.2054 
C = 6.8699 
D = .04604 

E = .02291 



By Bairstow's approximate factorization and one iteration for the phugoid 
roots, the quadratic is: 

X 2 + .00564X + .00332 = 0 

The period of the phugoid is 109 seconds and the time to damp to half 
amplitude is 246 seconds . 
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The short period roots were obtained from the remaining factor of the 
quartic : 

x 2 + 2 . 19976 X + 6.85419 = 0 

The period of the short period mode is 2.65 seconds and time to damp to half 
amplitude is O .63 seconds. 

The complete solution was obtained by using the Reeves Electronic Ana- 
log Computer and these results are shown in Fig. 3* The values obtained by 
the computer method and the stability quartic method are seen to be in excel- 
lent agreement and these values will be used as the basis for comparison of 
the various methods employed for simplification. 

SHORT PERIOD METHOD 

It is assumed that the short period mode is essentially a coupling be- 
tween the perturbations w and q alone. By this method, the equation is: 

AX 2 + BX + C = 0 

where : 

A = 1 
B = 2.201 
C = 6.84 

From this, the period is 2.66 seconds and the time to damp to half amplitude 
is O .63 seconds. This is in excellent agreement with the results obtained 
from the quartic. The computer solution is shown in Fig. 4. 

PHUGOID METHOD I 

The second method to be compared is the previously discussed phugoid 

simplification in which Aa was considered to be equal to JL . The quadratic 

Vn 
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equation for this method is: 



A\ 2 + B\ + C = 0 



where : 



A = 8.5306 
B = .01+336 

c = .02275 

The period by this method is 122 seconds and the time to one-half amplitude 
is 273 seconds. The period is seen to be about 12 percent in error. The 
computer solution is shown in Fig . 5 . 



PHUGOID METHOD II 



The third method to be compared is the method which requires that when 
Aa = 0, w = u tan a. The quadratic equation which results is: 

A\ 2 + B\ + C = 0 

where : 

A = 8.56 
B = .0555 
c = .031+ 

The method gives a period of 101 seconds and the time to damp to half ampli- 
tude is 213 seconds. The computer solution is shown in Fig. 6 for further 
comparison. It is seen that by this method the period is about eight percent 
too low. 



MODIFIED PHUGOID METHOD II 

The fourth method to be compared is the method evolved in this report. 
This method requires that when Aa = 0, w = u tan a. The quadratic equation 
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is: 



AX 2 + BX + C = 0 



where : 



A = 8,92 
B = .418 

c = .034 

The period by this method is 110 seconds and the time to damp to half ampli- 
tude is 29.6 seconds. The computer solution is shown in Fig. 7 for further 
comparison. It is seen that by this method the period is accurate to less 
than one percent . 



PHUGOID METHOD III 

The fifth and final numerical solution is the last phugoid case where 
w is considered to be zero and the second force equation is ignored. This 
method produces the cubic equation: 

AX 3 + BX 2 + CX + D = 0 

where : 

A = 1 

B = 1.60939 
c = .01788 

C = .01106 

This produces one negative real root and a complex pair. The period by this 
method' is 75. 7 seconds and the time to damp to half amplitude is 203 seconds. 
The computer solution is shown in Fig. 8 for further comparison. The error 
by this method is about 30 percent. 
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CASE II 



AIRPLANE DATA AND VELOCITIES 



Iyy = 80,000 SlUg ft 2 



V p = 279 ft/sec 
U x = 279 ft/sec 
W x = 31.9 ft/sec 
M = .25 



Q! - 6.5 degrees = 0 X 



S = 600 sq ft 
m = ll 80 slugs 



h = sea level 



X u = - 6.3 It sec/ft 
X v = 165 lb sec/ft 
= -162.6 lb sec/ft 
7 ^ = - 98 O lb sec/ft 
= - 13^5 lb sec 
= -102.7 lb sec 2 
1*^ = 153-7 lb sec 
Zq - -^325 lb sec 
Mq = - 115000 lb ft sec 



Computer solutions for the second case, Mach 0.25 at sea level are 
shown in Figs. 9 through 13 and the results are seen to parallel Case I. 

CONCLUSIONS 

Comparison of the numerical results and computer solutions of the 
several methods presented lead to the following conclusions: 

1. The concepts concerning the short period symmetric motion of an 
airplane when applied to the general equations provide an acceptable 
approximation to the period of that motion. 

2. Concepts concerning the period of the phugoid mode, when applied, 
provide accuracies ranging from 30 percent low to 12 percent high, 
depending upon the method of application. It was concluded, there- 
fore, that use of the existing simplifications does not accurately 
describe the period of the true motion because essential couplings 
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are neglected. 

3. The Reeves Electronic Analog Computer provided rapid accurate solu- 
tions to many sets of differential equations which had to “be solved. 
The utility of the computer in investigations such as this is unques- 
tionable. Variation of stability derivatives or coupling between 
equations could be accomplished by changing one or several resistances 
and the new solution was immediately obtainable . 
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